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ABSTRACT 

Context. Observations of the y-ray sky with Fermi led to significant advances towards understanding blazars, the most extreme class 
of Active Galactic Nuclei. A large fraction of the population detected by Fermi is formed by BL Lacertae (BL Lac) objects, whose 
sample has always suffered from a severe redshift incompleteness due to the quasi-featureless optical spectra. 

Aims. Our goal is to provide a significant increase of the number of confirmed high-redshift BL Lac objects contained in the 2 LAC 
Fermi/LAT cataloge. 

Methods. For 103 Fermi/LAT blazars, photometric redshifts using spectral energy distribution fitting have been obtained. The pho- 
tometry includes 13 broad-band filters from the far ultraviolet to the near-IR observed with SwiftfUVOI and the multi-channel imager 
GROND at the MPG/ESO 2.2m telescope. Data have been taken quasi-simultaneously and the remaining source-intrinsic variability 
has been corrected for. 

Results. We release the UV-to-near-IR 13-band photometry for all 103 sources and provide redshift constraints for 75 sources without 
previously known redshift. Out of those, eight have reliable photometric redshifts at z > 1.3, while for the other 67 sources we pro- 
vide upper limits. Six of the former eight are BL Lac objects, which quadruples the sample of confirmed high-redshift BL Lac. This 
includes three sources with redshifts higher than the previous record for BL Lac, including CRATES J0402-2615, with the best-fit 
solution at z a L9. 

Key words. Techniques: photometric, (Galaxies:) BL Lacertae objects: general. Galaxies: distances and redshifts 



^ ' 1. Introduction hole) (e.g., lAbdo etal.]|2010dl) to probes of the extra-galactic 

back ground li ght (EBL) through attenuation of y-ray photons 

d Smce Its launch m 2008, the Fermi Space Laboratory has dra- (e.g.. lAbdo et al. 2010c). One of the crucial parameters for these 

matically extended our view of the high-energy sky. The re- applications is the distance to the objects, which unfortunately is 
cendy released 24-month catalog (2LAC) of Active Galactic ^^^^ ^^^^^^^ ^^^^ ^^^^^ 
Nuclei (AGN) dete cted by the Large Area Telescope (LAT; 
lAtwood et al.l l2009l) revealed 885 high-sig nificance sources, 

the l arge majority of them being blazars dAckermann et alJ Two classes ofblazars dominate the 2LAC population. These 

l2011h . The latter form the most extreme class of AGN with are the flat-spectrum radio quasars (FSRQs, 310 sources) and 

their observational characteristics governed by the small an- BL Lac ob jects (395), named after the prototype BL Lacertae 

gle between their rela tivistic jets and the observer's sight (iHofi'meisteit 19291) . While for the former redshift measurements 

line dBlandford & Reesl ,1978). The resulting Doppler boosting are routinely performed using their strong emission lines at UV- 

makes blazars exceptionaUy bright sources at nearly all wave- optical wavelengths, the featureless, power-law opt ical spectra 

lengtiis and therefore visible out to high redshift. of BL Lac objects have proven to be a challenge (e.g.. [Shaw et all 

The scientific relevance of blazars is very broad, ranging 2009). Indeed, 220 of the 395 BL Lacs in the 2LAC (55 %) lack 

from laboratories for the physics and structure of relativistic jets redshift estimates. Until this incompleteness is r esolved, con- 

(and thus the extraction of energy from the central massive black elusions about the EBL, the blazar sequence (e.g., iFossati et al.l 
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et alJ[T998h . and the blazar population in gen- 
eral remain tentative at best. 

Several methods have been exploited to increase the BL Lac 
redshift sample. At low distance, one can utilize the remarkably 
uniform absolute brightness of t he giant elliptical BL Lac host 
galaxies dSbarufatti et all 120051; iMeisner & Romanil I20T0I) and 
very-high signal-to-noise optical spectroscopy to help to iden- 
tify weak emission or a bsorption features in a few other cases 
(e.g.. IShaw et al.l [2009h . An alternative method, applicable to 
the more distant sources, is the photometric redshift technique, 
which consists of fitting spectral energy distribution (SED) tem- 
plates to multi-band photometry. Neutral hydrogen along the line 
of sight to the blazar will imprint a clear attenuation signature at 
the Lyman limit and thus allows an accurate estimate of the red- 
shift of the absorber. Even though the absorber will be located 
somewhere along the line of sight and its redshift will thus not 
necessarily correpond to that of the blazar, photometric redshifts, 
Zphot, will provide a reliable lower limit on the blazar redshift. 

In this paper we explore the use of ultra-violet to near- 
infrared quasi-simultaneous photometry to measure the redshift 
of 103 2LAC blazars, 86 of them without previous redshift 
constraints. The sample selection is described in Sect. |2] while 
the observations and data reduction are detailed in Sect. [3] and 
Sect, m The results and their discussion are presented in Sect. |5] 
and Sect. |6] respectively. 



2. Sample Selection 

The sample is selected mainly from the clean sample 
of the 24-month catal og of AGN detected by Ferm//LAT 
(lAckermann et al.ll201 ih which includes 885 high-significance 
y-ray sources which are statistically associated with AGNs and 
located at high Galactic latitudes {\b\ > lOdeg). To accomodate 
the location of one of our follow-up instruments in the southern 
hemisphere, we applied a cut in declination (5j2ooo < +25 deg). 
Furthermore, the importance of the ultra-violet bands for the 
Zphot estimates lead us to discard all sources with Galactic fore- 
ground reddening o f Eb-v > 0.2 mag, as derived from the 
ISchlegel et al.l (Il998h maps. The final sample is composed of 80 
sources without known redshift but with o ptical or radio count er- 
part associations from Bayesian statistics (lAbdo et a i?'2010a b). 
Likelihood Ratio, or log N - log S methods (lAckermann et al.. 
l20Tlh and was further extended by eight sources that do not be- 
long to the clean 2LAC sampleQ In addition, observations for 16 
2LAC sources with existing redshift measurements are included 
for verification. We emphasize here that our sample is not sta- 
tistically complete, but biased by selection. The final list of 104 
targets is presented in Table [1] Of those, 82 have been classified 
as BL Lac objects, three as a FSRQ, and the remaining 19 are of 
unknown type (.Ackermann et aL,201 li) . 
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Fig. 1. Cumulative distribution of the (absolute) time between 
the mid-points of the GROND and ^wZ/f/UVOT observations. 
The median offset is 1.5 d. For longer UVOT observations, 
which were spread over several satellite orbits, the mid-time of 
the Mvw2-exposure was taken as reference. 



The UVOT pointings have been obtained as part of the Swift 
fill-in program from January 2010 to October 2011. With few 
exceptions, they comprise observations in three optical (m, h, v) 
and three UV (mvw2, uvm2, uvwl) lenticular filters, covering the 
wavelength range between 160 nm and 600 nm (IPoole et al.l 
2008). The exposure times vary from source to source de- 
pending on visibility and brightness of the optical counterpart. 
Approximate values are 100 s each in u, b, v and 200, 250, 400 s 
in uvwl, uvm2, uvwl, respectively. For a few sources multiple 
consecutive orbits, and thus longer exposures, were obtained. 

While UVOT needs to cycle through all filters in sequence, 
GROND observes in it's four optical g', r', i',z' and three near- 
IR J, H, Ks bands simultaneously. This capability is particularly 
important for studies of fast-varying sources such as blazars as it 
allows the construction of a reliable spectral energy distribution 
(SED) without the need to correct for source-intrinsic variabil- 
ity. Unfortunately, scheduling limitations rarely allowed simul- 
taneous observations with both instruments, UVOT from space 
and GROND from the ground. However, in order to minimize 
the impact of variability on the combined SED, GROND obser- 
vations were executed as close as possible in time to the Swift 
pointings. Except in a few cases where visibility or weather con- 
straints lead to offsets of more than lOd, observations were typ- 
ically achieved within 1-2 d of each other (see Fig.[T]). The Zphot 
uncertainties associated with these offsets will be discussed in 
Sect. 1431 

A typical GROND observation had an integration time of 
2.4 min in (g',r',i',z') and4.0min in {J,H, K^). 



3. Observations 

The observations were per formed with the U ltraviolet and 
Optical Telescop e (UVOT; 'Roming et al.l l2005h onboard the 
Swift satellite ( Gehrels et al. 2004~ and the Gamma-ray 
Optical/Near-Infrared Detector (GROND; lOreiner et alJ 120081) 
mounted on the MPI/ESO 2.2 m telescope at La Silla, Chile. 
Each source has been observed at least once with both instru- 
ments, although for some targets several epochs had been ac- 
quired. 



4. Data & Analysis 

4. 1 . Swift/UVOT and GROND 

UVOT photometry was carried out on pipeline processed sky 
images downloaded from th e Swift data center El following the 
standard UVOT procedure (iPoole et al.l l2008h . Source photo- 
metric measurements were extracted from the UVOT imaging 
data using the tool uvotmaghist (vl.l) with a circular source ex- 
traction region that ranged from 3'.'5-5" radius to maximise the 



' flagged due to analysis issues by the LAT team 



^ |http://www.swift.ac.uk/swift-portal| 
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Table 1. Target List (first five rows, rest in Journal only) 



Name 


FGL Name" 


Q'J200o'' 


OJ2000 


offset ["Y 


Eb-v [mag]'' 


CRATES JOOO 1-0746 


2FGLJ0000.9-0748 


00:01:18.00 


-07:46:26.9 


0.2 


0.03 


CRATES J0009+0628 


2FGLJ0009.0+0632 


00:09:03.93 


06:28:21.3 


1.1 


0.07 


CRATES J002 1-2550 


2FGLJ0021. 6-2551 


00:21:32.56 


-25:50:49.1 


0.8 


0.02 


IRXS J002209.2- 185333 


2FGLJ0022.2-1853 


00:22:09.26 


-18:53:34.9 


3.1 


0.03 


RX J0035.2+1515 


2FGLJ0035.2+1515 


00:35:14.71 


15:15:04.2 


3.4 


0.07 



": Fermi/LAT ID. 

*: GROND coordinates of the optical counterpart wi th typical uncertainties o f 0'.'3 in both directions. 
'^: Offse t of counterpart locati on from those given in lAclcermann et all d2011l) . 
'': From lSchlegel et al.l dl998h . 



signal-to-noise. In order to remain compatible with the effec- 
tive area calibrations, which are based on 5" aperture photome- 
try (Poole et al. 2008), an aperture correction was applied where 
necessary. This correction was at maximum 5-6 % of the flux, 
depending on the filter. For the further analysis, all magnitudes 
were converted into the AB system and are presented in Tab. |2l 
Typically achieved 3-cr limiting magnitudes are mvw2ab ~ 22.4 
and bAB ~ 19.6. 

The GROND data (Tab. O were reduced and analyzed 
with th e standard tools and methods described in iKriihler et alj 
(|2008|) . The g',r',i',z' photometry was obtained using point- 
spread-function (PSF) fitting and calibrated against observations 
of fie lds covered by the SDSS Data Release 8 ( Aihara et aTl 
1201 lb . Due to the undersampled PSF in the near-infrared, 
the J,H,Ks photometry was measured from apertures with 
the sizes corresponding to the Full- Width at Half Maximum 
(FW HM) of field stars and calibrated against selected 2MASS 
stars (ISkrutskie et al.ll2"006h . This resulted in Icr accuracies of 
0.04 mag (g' ,z'\ 0.03 mag (r', /'), 0.05 mag {J, H), and 0.07 mag 
(Ks) for the absolute calibration. While the SDSS calibration di- 
rectly provides magnitudes in the AB system, the near-IR pho- 
tometry required further transformation from the 2MASS -native 
Vega system into AB. Typically achieved 3-cr limiting magni- 
tudes are r^g ^ 23.5 and K^^ab ~ 19.8. 

Correction for Galactic foreground extinction was performed 
following the procedure describ ed in Cardelli et al. ( 1989) with 
Eb-v from lSchlegel et al.l (Il998l). For the UVOT ba nds, the cor- 



rection factors presented in 'Kataoka et al. ( 20081) were used 



Uncertainties in Eb-v (10 %; Schlegel et al .ll998h and in the red- 
dening law produce an additional contribution to the photometric 
error budget. As this systematic uncertainty is coupled between 
the the photometric bands for a given position in the sky, its im- 
pact on the SED fitting is smaller than the contribution to the 
photometry in each individual filter. An exact calculation is very 
complex. Thus, we adopt a conservative 5 % of the reddening 
value in each band which is added in quadrature to the photo- 
metric uncertainties. 



4.2. Counterpart selection and morphology 

Astrometric solutions for the GROND optical bands have been 
obtained through comparison with USNO or, if available, with 
SDSS measurements, achieving a typical rms of 0'.'3 in both co- 
ordinates. For the large majority of our targets (91) only a single 
optical source was found within 2" radius of the location given 
in the 2LAC catalog. For an additional twelve objects isolated 
candidate counterparts were detected at offsets of up to 6". In 
one case (CGRaBS 1407-4302), several blended sources made 
a reliable identification impossible. Counterpart coordinates and 




14 16 18 20 22 
r'-band Brightness [mag] 



Fig. 2. Observed GROND r'-band magnitude distribution of the 
optical counterparts. 



ofi^sets for the 103 sources with identification are given in Tab.[T] 
Most of the identified optical sources are bright and point-like 
with typical magnitudes in the GROND r'-band of 16-18 mag 
(Figure |2]l. 

The source morphology can be used as a prior for the se- 
lection of the SED template library (see Sect. 14.41 ) and redshift 
range, i.e., clearly extended sources will have a contribution 
from the host galaxy emission to the observed SED and are likely 
located at low redshift. A simple morphological distinction in 
point-like and extended has been performed on the GROND im- 
ages and will be used for discussing the reliability of the Zphot 
estimates. 

4.3. Variability correction 

The optical emission of blazars is known to vary on a wide 
range of time scales. Brightness variations range fro m a few 



tenths of a magnitude within minutes to days (e.g 



[T970I iMiller et al.1 [19891 ICarini et al.l [1991 lUrrv et al 



Racind 
I993h 



to several magnitudes over weeks to years (e.g.. ' Ciprini et al.) 
2003; Raiteri et al. 2005). This variability contributes signifi- 
cantly to the uncertainties when constructing an SED from non- 
simultaneous multi-band, multi-instrument observations. 

As GROND observes in all seven bands at the same time, the 
g' to Ks photometry can be considered as a snapshot of the SED 
and is thus unaffected by variability on time scales longer than 
the exposure of an individual observation. This leaves two ar- 
eas where a proper treatment of variability has to be performed. 
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Table 2. 5w;/f/UVOT photometry (first five rows, rest in Journal only) 



Name 


UT Date" 








AB Ma 


gnitude* 
















i(vw2 


uvm2 


uvwl 


u 




b 


V 




[mag] 


CRATES JOOOl... 


2010/10/14 10:40 


18.99 ± 0.08 


18.74 ±0.10 


18.58 


±0.09 


18.18 = 


t0.09 


17.88 ±0.10 


17.90 ± 


0.19 


0.01 


CRATES J0009... 


2010/12/07 08:55 


20.16 ±0.15 


19.71 ±0.19 


19.69 


±0.18 


19.24 = 


t0.19 


18.51 ±0.18 


18.14 ± 


0.24 


-0.24 


CRATES J0021... 


2010/11/18 18:32 


18.86 ±0.08 


18.90 ±0.27 


18.57 


_i_ n no 

± u.uy 


18.15 = 


tO.lO 


17.82 ±0.11 


17.93 ± 


0.21 


0.09 


IRXS J002209.2... 


2010/11/18 13:49 


18.36 ±0.07 


17.95 ±0.08 


17.91 


-t- n OS 

± U.Uo 


17.48 = 


t0.07 


17.20 ±0.08 


16.97 ± 


0.12 


-0.08 


RX J0035.2... 


2011/01/14 01:08 


18.36 ±0.08 


17.95 ±0.08 


17.91 


± 0.08 


17.48 = 


t0.07 


17.20 ±0.08 


16.97 ± 


0.12 


-0.24 


Approximate start time of uvwl exposure. 






















: Corrrected for Galactic foreground reddening. Upper limits are 3-cr. 


















Variability-correction factor to be applied to GROND photometry (see text). 
















Table 3. GROND photometry (first five 


rows, rest in Journal only) 


















Name 


UT Date" 










AB Mag 


— '■ 1 — TT" 

nitude 














8' 


r' 






z 




J 


H 






CRATES JOOOl... 


2010/10/14 03:33 


17.80 ± 0.05 


17.43 ± 0.05 


17.18 


±0.05 


16.86 = 


t0.05 


16.48 ±0.05 


16.16 ± 


0.06 


15.76 ± 0.08 


CRATES J0009... 


2010/12/12 00:41 


18.91 ±0.05 


18.49 ±0.05 


18.15 


±0.05 


17.89 = 


t0.06 


17.34 ± 0.06 


16.90 ± 


0.06 


16.46 ± 0.08 


CRATES J0021... 


2010/11/23 00:13 


17.69 ± 0.05 


17.43 ± 0.05 


17.20 


±0.05 


16.98 = 


t0.05 


16.76 ±0.05 


16.51 ± 


0.05 


16.15 ±0.07 


IRXS J002209.2... 


2010/11/23 00:21 


17.30 ±0.05 


16.92 ±0.05 


16.71 


±0.05 


16.51 = 


t0.05 


16.22 ±0.05 


16.00 ± 


0.06 


15.57 ± 0.08 


RX J0035.2... 


2011/01/14 00:48 


17.39 ±0.05 


17.10 ±0.05 


16.96 


±0.05 


16.77 = 


t0.05 


16.41 ±0.06 


16.17 ± 


0.05 


15.92 ± 0.08 



": Exposure start time. 

*: Corrrected for Galactic foreground reddening. Not corrected for variability. Upper limits are 3-it. 



i) between the exposures in the individual UVOT filters and ii) 
between the GROND and UVOT pointings. 

UVOT observes in each band separately and, for our pro- 
gram, typically cycled through its filters in a specific order 
(uvwl , u, b, uvwl, V, uvmT). A complete sequence in all six bands 
took ^ 12min, during which time the target may vary in bright- 
ness. Without a priori knowledge of the spectral shape, this vari- 
ability can only be accounted for statistically, i.e., by including 
an additional contribution to the systematic uncertainties. For 
this purpose we analyzed the photometry for those blazars for 
which multiple UVOT exposures in a single band were available. 
Here, we focused on the bluest filter {uvw2) which would be 
least affected by a potential host galaxy contribution and should 
thus give a good representation of the maximum variability. In 
Fig.[3]the distribution of variability as function of time between 
two uvwl exposures is shown for 25 targets with accurate pho- 
tometry (dniunvi < 0. 1 mag) in each image. While variability as 
large as Aniuvwi ~ 0.4 mag can be seen for some sources also on 
time scales of Af < 0.1 d, the median variation for Af < 50d is 
K 0.1 mag, comparable to the photometric uncertainty. 

The shortest time sampled between two uvwl exposures in 
this sample is one satellite orbit (x 90min) and therefore longer 
than it takes to cycle through the six filters. However, it is rea- 
sonable to assume that the median variation on a 12min time 
scale is not much larger than the one over 96 min. Thus, we con- 
servatively adopted Am = 0.1 mag as an additional contribution 
to the systematic uncertainties for all UVOT filters. 

In order to correct between GROND and UVOT we made 
use of the spectral overlap provided by the two instruments 
(lKruhIeretalJl20Ilh . Under the assumption that the SED re- 
mains unchanged and can be approximated by a power law and 
thus, that the GROND and b, v photometry follow the same spec- 
tral slope Q, color-terms can be used to derive the normaliza- 



^ This is reasonable for redshift z < 3.5 and for sources where the 
host galaxy contribution to the emission is negligible. See Sect. l4.4l for 
further discussion. 




1.0 10.0 
At [days] 



Fig. 3. Absolute brightness changes in the uvwl filter as func- 
tion of time between two UVOT pointings for all 25 sources 
with multiple observations. The black crosses show the median 
values (and median absolute deviation) and the grey dots mark 
all 174 pairs of observations. While micro-variability (Am„vu 2 ~ 
0. 1 mag) dominates on time scales shorter than 50 d, signif- 
icantly larger variability (Ato„,,,j,2 > 0.3 mag) is observed at 
At > 50 d. 



tion offset. Here we used the calibration with the largest spectral 
overlap, namely the one between b, g', and r': 



b-g' ^O.n{g'-r') + 0.03{g'-r'f 



(1) 



which is vahd for -1 < (g' - r') < 1 (iKriihler et al .11201 ih . 
These corrections, AmagcR^uv (Tab |2]i, were calculated for all 
sources with accurate photometry in b, g', and r' (6m < 0.2 mag). 
They were computed for each object individually and applied to 
the GROND photometry, shifting it in line with UVOT. Typically 
values are of the order of +0.1 mag (Fig. |4]i. In cases where 
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Fig. 4. Absolute brightness offset as function of time between the 
GROND and Swift/UYOT observations. The black crosses show 
the median values (and median absolute deviation) and the grey 
dots mark all 90 sources for which the offset was determined 
(excluding eleven sources with large photometric uncertainties, 
missing b, g', or r' photometry, or strong host galaxy contribu- 
tion). 



multiple epochs were available, we used the closest pair of 
GROND and UVOT observations. The final SED for each source 
corresponds to the best possible reconstruction between 160- 
2200 nm at the time of the UVOT observation. No correction 
was performed for sources where the photometric uncertainty 
(6m > 0.2 mag in b,g', and/or r') was larger than the typical 
AmagGR^vv- 

For a few sources, the emission from the underlying 
host galaxy contributes significantly to the overall SED (see 
Sect. 14.41 ). While the blazar emission can still vary, the host con- 
tribution naturally remains unchanged, violating the assumption 
of a constant SED slope. Thus, no variability correction was per- 
formed in those cases where the UV to optical SED clearly de- 
viates from a power law. 

4.4. SED fitting 

The photometric redshift computation wa s performed with the 
pubhcally ava ilable LePhare code v.2.2 jArnouts et alj [19991 
lllbert et al.l20"0 6). The program uses the simple fitting method 
to differentiate between a theoretical and an observed photo- 
metric catalogue. We used three custom template libraries for 
comparison with our GROND and UVOT photometry. The first 
was composed of 40 power law SEDs of the form oc A'^ 
with /3 ranging from to 2 in steps of 0.05 and describes pri- 
marily the central-engine dominated blazars. Here, the source- 
intrinsic extinction was assumed to be zero and the luminos- 
ity prior was choosen to cover th e expected range of absolut e 
magnitudes (-20 < Mg, < -30: IVeron-Cettv & Veronll2010h . 
Such a simple model is a good approximation for the UV-near- 
IR SEDs of BL Lacs, which form the bulk of our sample. Also 
FSRQs, if their emission lines are not dominant, can be mod- 
eled with a power law to first order The second library con- 
tains te mplates of normal (inactiv e) galaxies and galaxy/AGN 
hybrids ( ISalvato et al.ll2009l 1201 I h and has been implemented in 
order to model host galaxy dominated objects. If the host con- 
tribution is significant, the 4000 A break can emerge and poten- 
tially be misinterpreted as a Lyman-limit. While for the galaxy 



models intrinsic extinction was again neglected, a different lu- 
minosity prior of -8 < Mg' < -30 was chosen to allow for the 
expected fainter host magnitudes. The third libr ary consists of a 
wide range of stellar templates dPickle s 1998: B ohlin et al.ll995l : 
IChabrier et alJi2000,) . included to test against potential false as- 
sociations. For each source in our sample, all three libraries were 
fit independently. 

The most prominent spectral feature used to measure Zphot 
for blazars is caused by the absorption through neutral hydrogen 
along the line of sight. At z 0.8 the Lyman-limit moves into 
the wavelength covered by the UVOT and starts to suppress the 
flux in the uvw2 filter While this is the theoretical lower limit on 
Zphot for the combined GROND and UVOT coverage, uncertain- 
ties associated with the photometry and modeling will shift it to 
a higher redshift. Usually one would use a spectroscopic train- 
ing sample to evaluate the redshift range and accuracy which is 
accessible with our SEDs. As the majority of the sources with 
known redshift in our sample are located at z < 1, they alone 
do not suffice to provide a statistically meaningful test of the 
method. Instead we assess the reliability in our photometric red- 
shift computation with a Monte-Carlo approach. 

Here, we simulated 27,000 test power-law SEDs with spec- 
tral slopes ranging from 0.5 to 2, redshifts from to 4, and at 
three apparent magnitudes, r' = 17, 18, 19 mag. For each source, 
the photometry in the individual bands was scattered around the 
model magnitudes with a brightness-dependent statistical con- 
tribution and with a component representing the systematic un- 
certainties from the photometric calibration and variability cor- 
rection. The resulting SEDs were fed back into LePhare and the 
output Zphot was compared to the input value. 

As can be seen in Fig.|5] there is a good correspondence be- 
tween the input and the recovered redshift for Zsim > 1-2. In par- 
ticular for bright sources (r' < 18 mag) the number of outliers t], 
defined as the number of sources with |Az(l -i-Zsim)l > 0.15, drops 
steeply above this redshift. For fainter sources (r' w 19 mag), the 
increasingly less constrained UV photometry causes the outlier 
fraction to remain high (w 20 %) out to Zsim ~ 2. A more quan- 
titative selection can be made when using P^, the integral of the 
probability distribution function J f(z)dz at Zphot + 0.1(1-1- Zphot), 
which describes the probability that the redshift of a source is 
within 0.1(1 + z) of the best fit value. When choosing a cut at 
Pz > 90 %, the majority of the outliers at Zphot > 1-2, as well as 
nearly all the solutions with Zphot < 1-2 disappear (inset Fig.|5]l. 
This confims that for power-law sources with z < 1 .2 the avail- 
able photometry can only place an upper limit, and suggests that 
Pz > 90 % can be used as a criterion to identify reliable photo- 
metric redshift solutions. Selecting a higher threshold, e.g., 
99 %, can further reduce the outlier fraction but will also shift 
the lower constrainable redshift bound to higher values. In order 
to reach such a very tightly peaked redshift probability distribu- 
tion function more than one photometric band has to be affected 
by the Lyman-limit, limiting the redshift range to z ^ 1-5. 

5. Results 

The fit results for all 103 sources with identified counterparts 
are presented in Tab. |4] Here we give the photometric redshifts 
(and their 90 % confidence errors), the Pz values, and the best 
fit models for the power law and galaxy templates. No source 
requires a stellar template. 

As discussed in the previous section, Pz > 90 % can be used 
as a good reliability criterion for the photometric redshift solu- 
tion. Applying this cut to the power law model fits results in 
a sample of 15 sources with < 30. Next, we also require 
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Table 4. Excerpt of the results table. Only entries for those sources with reliable redshifts of Zphot > 1 -2 are shown, (rest in Journal 
only) 



Name 



Sphot,best 
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-om 
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model 


8.7 


75.5 


12249 1.60_TQSO 1.40 


13.6 


80.5 


pLQSOJDR2_029_tO 


11.2 


38.4 


12249 1.50_TQSO 1.50 


13.0 


62.3 


12249 1.60.TQSO 1.40 


9.5 


48.9 


12249 1.60.TQSO 1.40 


21.5 


99.6 


pl.l22491.30.TQSO1.70 


7.5 


58.1 


12249 1.60.TQSO 1.40 


24.6 


83.7 


12249 1.60.TQSO 1.40 


11.9 


94.9 


pl.l22491.30.TQSO1.70 


18.4 


100.0 


12249 1.60.TQSO 1.40 


7.1 


64.2 


pl.l22491.20.TQSO1.80 



RX J0035.2+1515 1 

PKS 0047+023 1 

PKS 0055-328 1 

PKS 0332-403 1 

CRATES J0402-26 15 1 
SUMSSJ053748 571828 1 

PKS 0600-749 1 

CRATES J0630-2406 1 

OM 235 1 

CRATES J1312-2156 1 

CLASS J2352+ 1749 1 



1.426 



1.549 
1.491 



Best photometric redshift, see text. 

Spectroscopic or imaging redshift (if known) from lAckermarm et alj ( 120 iTh . 
^: Photometric redshift with 90 % confidence uncertainties. 
'': Redshift probability density at Zpi,ot ± 0.1(1 + Zphot)- 
'': Spectral slope for power law model of the form A^^. 

^: Starburst/QSO hybrid templates with varying contributions of the two components. Template names reflect the relative contribution of the 
starburst (122491.NN) and AGN (TQS01.NN) models to the hybrid, "pi" marks an additional power-law component at short wavelengths. See 
[Salvato et al., c201 1.) for details on the templates. 



that a source can not be fit well with a galaxy/AGN-hybrid tem- 
plate at low redshift (z < 1.2). This step is required in order 
to identify those sources for which a degeneracy in the pho- 
tometric redshift solution does not allow a distinction between 
a low-redshift galaxy/AGN-hybrid and a blazar at higher red- 
shift. In analogy with our selection for a good redshift constraint 
(Pz > 90 %), we consider Pz,gai < 90 % as a reliable criterion 
that no adequate low-redshift galaxy/AGN-hybrid solution was 
found. This removes further seven sources from the sample and 
leaves eight candidates, all unresolved in the GROND images 
and with Zphot,pi > 1.2. 

For three additional sources the power law and galaxy mod- 
els fit similarly well (Pz.pi ~ Pz.gai > 90 %) and give compara- 
ble high-redshift solutions, indicating that the Lyman-limit was 
fit by power law and galaxy templates alike. Including those, 
we end up with a sample of eleven sources with reliable photo- 
metric redshift computation, all , with Zphot ranging from Zphot ~ 
1.28 to Zphot ~ 1 -92. This encompasses eight sources without 
previous redshift measurements and three blazars with known 
redshift. For the latter, the photometric redshifts are within 3- 
cr of the known value (see Fig. |5]l, indicating an accuracy of 
Az/(1 + Zspec) < 0.12 for our Zphot method described in Sect. 14.41 
The fit results and SEDs for the eleven sources are presented in 
Tab.|4]and Fig.|6] respectively. 

For 21 sources, the galaxy library provides the better fit with 
Pz.gai > 90 %, Pz.pi < 90 %, and ^l^i < 30. All of these are 
best described by starburst/QSO hybrid templates at z < 1.3 
and none require an elliptical galaxy, i.e., a strong 4000 A break. 
The most significant difference between the starburst/QSO hy- 
brid templates and the power-law models is that the former ex- 
hibit strong e mission lines originating from star-formation (see 
ISalvatoetaP 120091 for more details on the templates). These 
lines can be fit to small deviations in the photometry from a 
power-law and thus provide low values and apparently well- 
constrained Zphot solutions. However, as described in Sect. 14.31 
the variability correction is not reliable when applied to non- 
power-law sources, as the color terms, and thus the correction 
factors, will be erroneous. Also, 18 of these sources have pre- 



viously been classified as BL Lac objects and are therefore 
unlikely to show strong emission lines in their optical spec- 
tra. Given these considerations, we do not regard their galaxy- 
template photometric redshifts to be reliable. 

There are 79 sources which are well-fit by power-law tem- 
plates (x^ < 30) but have a low Pz. In other words, the spec- 
tral slope is well defined but the photometric redshift cannot 
be unambiguously determined. As expected, the best solutions 
are predominantly at low redshift (z < 1.2) with typical red- 
shift probability distributions that are flat down to z = 0. For 
those we give the 90 % upper limit of the photometric redshift.. 
This sub-sample also includes eleven sources with a previously 
known redshift from spectroscopy or imaging, all in agreement 
with our limits. 

We also provide the 90 % upper limit for four sources for 
which the power-law and galaxy libraries give similarly good 
constraints (Pz.pi ~Pz,gai > 90 %) but with very different red- 
shift solutions. This degeneracy can occur, e.g., when a break 
in the SED is similarly well fit by the Lyman limit of a power- 
law model and with the 4000 A break of a galaxy template. No 
Zphot.best is givcn for ten sources for which no satisfactory fit 
(X^ > 30) was obtained. 



6. Discussion & Conclusion 

In this paper we presented redshift constraints for 103 blazars 
from the 2LAC catalog using UV-to-near-IR multi-band pho- 
tometry obtained quasi-simultaneously with SwiftfUWOT and 
GROND. We provided the first reliable redshift measurements 
for eight sources and new upper limits for an additional 66 tar- 
gets. Of the eight sources with reliable redshift, seven are lo- 
cated at Zphot > 1-3. Six of those belong to the BL Lac popula- 
tion. For comparison, out of the total 395 BL Lac in the 2LAC 
sample only two sources have previously been known to lie at 
z > 1.3 (Fig.llli. Redshifts for these two, PKS 0332-403 and 
CRATES J13 12-2156, have also been confirmed with our pho- 
tometric observations. 
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Fig. 5. Recovered best-fit photometric vs input redshift (top), ac- 
curacy (middle), and outlier fraction (bottom) for 27,000 sim- 
ulated sources with spectral slope of 0.05 < /3 < 2 and op- 
tical brightness of r' = 17, 18, 19mag. The pink circles indi- 
cate the location of three blazars in o ur sample with know n 
redshift, PKS 0332-403 (Zspec = 1.351. iBergeron et alllloTlh . 
OM 2 35 (z.np. = 1.549), and CRATES Jl 3 12-2 156 (zspec = 
1 .49 1 . lAckermann et all I20TTI). as well as that of PKS 0537- 
286 rz.p.. = 3.104. IWrightetanil978h . a LAT-non-detected 
blazar with available GROND and UVOT photometry. The green 
area in the middle panel corresponds to the forbidden parameter 
space where Zphot > 0. Systematics are visible at Zphot < 0.8. 
The large majority of the fits have < 30 (for typically 10 
d.o.f.) with the exceptions shown in light grey. The inset in the 
top panel displays the results if only solutions with > 90 % 
are selected. 



The six new BL Lac redshifts at z > 1.3 represent a dra- 
matic (from two to eight) increase of the confirmed high-z Fermi 
sample and thus demonstrate the opportunity that the SED tem- 
plate fitting technique holds for obtaining photometric redshifts 
for BL Lac sources. For our sample, this was possible due to 
the densely-covered wide spectral range (1 60-2200 nm), neces- 
sary for a reliable constraint of the spectral slope, the excellent 
ultra-violet coverage to measure the Lyman-limit, and the quasi- 
simultaneity of the observations, important for minimizing the 
impact of source-intrinsic variability. The method applied in this 
work overcomes significant challenges inherent in other redshift 
techniques, namely the simplicity of the optical emission of BL 
Lacs manifested as the power-law-shaped synchrotron spectral 
component. This complicates spectroscopic redshift measure- 
ments, which are in most cases limited to the optical wavelength 
regime and thus insensitive to the Lyman-limit. Instead, spec- 
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Fig. 6. SEDs for the eleven sources for which photometric red- 
shifts larger than 1 .2 were estimated (see Tab. |4]i and one ex- 
ample for Zphot < 1-2. From bottom to top: CRATES J0402- 
2615, CRATES J1312-2156, OM 235, CRATES J0630- 
2406, SUMSS J053748-571828, PKS 0600-749, PKS 0332- 
403, CLASS J2352+1749, PKS 0047+023, PKS 0055-328, 
RX J0035.2+1515, and BZB J0543-5532. 



troscopy relies on the detection of very faint emission features, 
mainly from the underlying host galaxy. It therefore requires 
very high signal-to-noise, often at the cost of long exposure 
times. On the other hand, the photometric redshift method for 
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Fig. 7. Redshift distribution for the nine BL Lac objects with re- 
liable photometric redshifts (blue filled histogram) together with 
the distri bution for the BL Lac in the 2LAC catalog (empty his- 
togram; Ackermann et alj201 lb . Sources with z < 0.5 have been 
omitted for clarity. 
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BL Lac does not require significant telescope time, in particular 
when data from efficient multi-channel instruments like GROND 
are available. Its weakness, however, is that the method relies 
on the detection of a particular spectral feature, which, albeit 
strong, is located in the rest-frame far-ultra-violet. Forz < 3, this 
is only accessible with space-based ultra-violet telescopes, and 
even then usually limited to redshifts above z ~ 1.2, as shown 
above. 

The sensitivity of this method to high-redshift sources also 
allows the placement of upper limits for those SEDs that do not 
show the imprint of the Lyman-limit. Only for three targets is 
this boundary above z = 2. In one case, CRATES J0250+1708 
(z < 3.1), the optical counterpart is too faint (r' » 20.7 mag) 
and only upper limits in all six UVOT bands could be obtained. 
The counterpart of CRATES J0705-4847 is only detected in the 
GROND J (w 19.7 magAB) and H bands and its SED and redshift 
are thus poorly constrained. Finally, the SED of ATG20 JO 124- 
0625 (z < 2.46) shows a significant break at z ~ 1.9 + 0.5. 
However, as the redshift probability distribution is broad, and 
the of the power law fit comparable to that of a z ~ 0.4 so- 
lution for a galaxy template, the photometric redshift is consid- 
ered to be unreliable. Except for those three, no other source in 
our sample has a best-fit photometric redshift at z > 2. Thus, 
CRATES J0402-2615 can now be considered the most distant 
known BL Lac with a measured redshift of z ~ 1 .92. 

Fig.|7]indicates that our high-redshift findings are the natural 
extension of the existing 2LAC redshift sample. However, due to 
the incompletenesses of both samples, we refrain from drawing 
any quantitative conclusions at this stage. We note, however, the 
our r esult is genera l agree ment with the theoretical predictions 
from lGiommi et al.l (1201 11) . A more detailed physical interpreta- 
tion of the establishment of an increased fraction of high-z BL 
Lacs in the 2LAC sample will be reported separately. 
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